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摘 要 :近年 来 ,温度 植被 干旱 指数 (temperature vegetation drought index , TVDI) 在 遥感 旱情 监测 中 得 到 广泛 的 应 用 
但 传统 的 TVDI 没有 考虑 地 形 因 子 对 指数 精度 的 影响 。 本 文 以 陕西 省 为 例 , 利 


J DEM 数据 对 地 表 温 度 做 高 程 修 


正 ,同时 引入 归 一 化 差 值 山 地 植被 指数 (normalized difference mountain vegetation index,NDMVI) 代替 归 一 化 植被 指 
数 (normalized difference vegetation index,NDVI) ,以 减弱 地 形 起 伏 对 植被 指数 的 影响 ,进而 构建 新 的 TVDI。 分 别 利 
用 修正 前 后 的 地 表 温 度 以 及 NDMVI 和 NDVI 构建 了 4 种 TVDI, 并 分 析 4 种 指数 与 土壤 湿度 及 降水 的 关系 ,分 析 比 
较 其 监测 农业 干旱 的 精度 。 结 果 表 明 : 高 程 与 地 形 起 伏 对 TVDI 监测 干旱 的 精度 有 较 大 的 影响 。 经 过 修正 后 的 
TVDI 与 传统 TVDI 相 比 ,其 与 土壤 相对 湿度 的 相关 系数 从 -0. 342 提高 至 -0.711 ,上 且 能 更 准确 地 体现 降水 的 变化 
特征 ,说 明 修 正 后 的 TVDI 能 够 更 有 效 地 反映 区 域 土壤 水 分 状况 ,可 对 农业 干旱 进行 更 加 精确 地 监测 。 


关键 词 :干旱 监测 ; 地 形 因 子 ; TVDI; NDMVI; 修正 ; 陕西 


近年 来 , 随 着 气 修 变 暧 等 全 球 性 环境 问题 的 凸 
显 , 极 端 灾害 事件 频繁 发 生 '""。 干旱 ,尤其 是 重大 


干旱 灾害 已 经 威胁 到 社会 和 国家 稳定 ,干旱 的 发 生 
与 持续 对 农业 生产 .粮食 安全 生态 环境 以 及 国民 经 
济 等 产生 严重 的 影响 ”) 。 因 此 ,关于 干旱 ,尤其 是 
农业 干旱 的 监测 研究 受到 了 学 者 们 的 广泛 关注 ” 。 
传统 的 农业 旱情 监测 主要 是 通过 土壤 墙 情 站 点 或 气 
象 站 点 ,获取 农田 土壤 湿度 和 降水 量 等 信息 来 描述 
旱情 状况 , 虽 能 真实 地 反映 干旱 发 生 状况 ,但 监测 精 
度 取决 于 地 面 站 点 的 分 布 密 度 , 较 难 详细 反映 农业 
干旱 状况 , 且 览 测 成 本 较 高 "“”。 而 遥感 技术 能 快 
速 获 取 地 表 和 覆盖 、 辐 射 条 件 等 信息 ,具有 观测 范围 
广 \ 获 取信 息 量 大 、 速 度 快 . 实 时 性 好 动态 性 强 等 优 
点 ,已 成 为 目前 农业 旱情 监测 普遍 采用 的 方法 "。 
其 中 温度 植被 干旱 指数 (TVDI) 因 模 型 建立 简单 、 监 
测 精度 较 高 ,在 遥感 干旱 监测 中 得 到 了 广泛 的 研究 
和 应 用 。 

TVDI 是 基于 光学 与 热 红外 遥感 通道 数据 进行 
植被 覆盖 区 域 表 层 土壤 水 分 反 演 的 方法 来 进行 干旱 
监测 的 ,以 卫星 得 到 的 地 表 温度 T, 数据 与 NDVI 为 
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基础 进行 计算 ,然而 地 表 温 度 与 NDVI 分 别 受到 高 
程 和 地 形 起 伏 的 影响 。 一 方面 , 随 着 高 程 的 增加 大 
气温 度 逐 渐 降 低 , 气 温 通 过 满 流 的 形式 与 地 表 温 度 
进行 热 交 换 , 从 而 对 地 表 温 度 产 生 影响 。Ran 等 
用 高 程 对 TVDI 模型 的 地 表 温 度 进行 了 修正 , 结 

表明 经 过 高 程 修正 后 获取 的 TVDI 值 与 实测 土壤 相 
对 湿度 具有 更 好 的 相关 性 。 刘 立 文 等 在 研究 多 
种 植被 指数 与 地 表 温 度 的 基础 上 ,引入 DEM 数据 对 
地 表 温度 做 地 形 修 正 , 得 出 地 表 温 度 在 经 过 地 形 修 
正 后 能 够 更 好 地 反映 研究 区 的 地 表 实 际 状况 。 另 一 
方面 ,遥感 植被 指数 是 通过 多 光谱 遥感 数据 波段 的 
线性 或 非 线 性 组 合 来 描述 地 表 植 被 状况 的 一 种 光谱 
特征 度量 指标 心 -”] 。 在 植被 覆盖 区 域 ,地 形 起 伏 会 
导致 愧 感 影像 中 同 种 或 相似 地 物 呈 现 不 同 的 光谱 反 
射 特征 ,这 种 反射 特征 的 差异 也 会 给 植被 信息 估算 
带 来 误差 ,使 得 分 布 在 地 形 区 域 较为 复杂 的 植被 信 
息 很 难 完全 提取 "”。 地 形 起 伏 引 起 的 光谱 信息 变 
化 严重 影响 了 植被 信息 的 提取 精度 ,成 为 植被 遥感 
害 息 准确 解 译 的 主要 障碍 '"“""“ 。 为 消除 或 者 抑制 
地 形 的 影响 , 吴 志 杰 等 "利用 山地 植被 指数 提取 植 
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被 信息 ,并 取得 良好 的 效果 。 综 上 ,高 程 及 地 形 起 估 
分 别 对 人 和 NVDI 造成 了 较 大 的 影响 ,然而 , 以往 
在 TVDI 干旱 监测 研究 中 并 没有 将 这 种 影响 完全 考 
虑 进去 。 

为 了 消除 和 减弱 高 程 和 地 形 起 伏 对 TVDI 计算 
结果 的 影响 ,提高 其 监测 农业 旱情 的 精度 ,本 文选 择 
地 形 起 伏 较 大 的 陕西 省 作为 研究 区 ,在 分 析 TVDI 
特征 空间 的 基础 上 ,引入 DEM 数据 对 7, 进行 高 程 
修正 ,同时 引入 NDMVI 来 代替 NDVI 反映 植被 信 
息 。 通 过 修正 后 的 地 表 温 度 和 NDMVI 构建 新 的 温 
度 植被 干旱 指数 ,利用 农田 土壤 相对 湿度 数据 及 降 
水 数据 评价 不 同 指数 监测 旱情 的 效果 ,以 期 提高 温 
度 植被 干旱 指数 (TVDI) 反 演 土 壤 水 分 的 准确 性 ,为 
遥感 干旱 监测 提供 新 的 思路 与 方法 。 


1 资料 与 方法 


1.1 研究 区 概况 
陕西 省 位 于 我 国 西北 内 陆 , 地理 位 置 为 
105°29' ~111°15'E 31°42’ ~ 39°35'N ,地 域 南北 长 ， 
RTF , 横 跨 黄河 和 长 江 两 大 流域 中 部 (图 1)。 
BETE Soe AL LL BEA ,河川 纵横 ,地 势 总 体 特点 
为 南北 高 ,中 部 低 ,海拔 分 布 落差 较 大 ,高程 范围 为 
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1 研究 区 概况 示意 图 
Fig. 1 Geographical location of the study area 
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171 ~3 621 m, 以 北山 与 秦岭 为 界 , 全 省 自 北向 南 分 
为 陕 北 高 原 、 关 中 平原 和 秦 巴 山地 三 大 地 貌 区 "。 
陕西 省 境内 气候 差异 很 大 ,由 北向 南 渐次 过 渡 为 温 
带 \ 暖 温带 和 亚热带 ,整体 属 大 陆 季 风 性 气候 ,因而 
四 季 分 明 , 雨 热 同期 , 易 发 生 干 旱 '“ 。 全 省 降水 量 
分 布 不 均 , 自 东南 向 西北 地 区 递减 ,年 平均 降水 量 为 
300 ~1 000 mm"? 。 陕 北 地 区 气候 较为 干燥 ,干旱 
发 生 频 率 高 达 73. 3% OP ;关中 地 区 为 半 湿 润 气 候 
IX, FRAC PR ee OM ; 陕 南 地 区 旱灾 时 有 发 生 , 其 中 
商洛 地 区 易 受 旱灾 影响 。 根 据 省 气象 局 相关 资 
料 记载 显示 ,2010 年 陕西 省 出 现 较为 严重 的 春 夏 旱 
情 , 对 农作物 的 生长 产生 严重 影响 ,其 中 旱情 较 重 的 
渭南 .成 阳 两 地 作物 受 旱 面积 已 超过 2. 0 x 10° 
hm? 23) i 
1.2 数据 介绍 

本 研究 采用 的 遥感 数据 来 源 于 NASA 的 MO- 
DIS 产品 数据 (MOD11A2 和 MOD13A2 ,时 相 是 2010 
年 4 月 7 日 至 9 月 30 日 ), 行 列 号 为 h26v05 与 
h27v05。MOD11A2 是 8 d 合成 的 地 表 温 度 产品 ， 
MOD13A2 是 16 d 合成 的 植被 指数 产品 ,在 获得 人 研 
究 区 域 地 表 温 度数 据 与 植被 指数 数据 之 后 ,将 
MOD11A2 数据 采用 均值 法 合成 为 16 d 的 地 表 温 度 
数据 。DEM 采用 30 m 的 ASTER GDEM 数据 , 重 采 
样 为 1 km 分 辩 率 的 数据 。 

降水 数据 来 源 于 中 国 地 面 气候 资料 日 值 数据 集 
中 的 每 日 20 时 至 次 日 20 时 降水 记录 资料 ,2010 年 
陕西 省 共有 20 个 降水 记录 站 点 ,其 中 有 效 站 点 为 
17 个 。 土 壤 墙 情 数 据 来 源 于 中 国 农作物 生长 发 育 
和 农田 土壤 湿度 旬 值 数据 集 , 研 究 采 用 该 数据 集中 
陕西 省 28 个 农业 气象 站 点 测 得 的 20 cm 土壤 相对 
湿度 , 当 土 壤 水 分 饱和 时 土壤 相对 湿度 记 为 99% ， 
主要 为 降水 或 者 人 为 灌溉 导致 ,为 了 保证 使 用 的 精 
度 ,将 土壤 相对 湿度 为 99% 的 值 剔 除 。 两 个 数据 集 
均 下 载 自 中 国 气 象 数据 网 (http : //data. ema. en) 。 
1.3 研究 方法 
1.3.1 TVDI 和 干旱 监测 原理 ”Sandholt 等 * 在 研究 
土壤 湿度 时 发 现 , 如 果 研 究 区 的 植被 覆盖 包含 从 裸 
地 到 植被 全 覆盖 ,土壤 湿度 包含 从 极 干 旱 到 极 湿润 
的 各 种 情况 ,那么 ,以 遥感 获得 的 NDVI 和 地 表 温 度 
CT.) 为 横 纵 坐标 组 成 的 散 点 图 呈 三 角形 的 特征 (图 
2), H NDVI ~T, 特征 空间 中 存在 许多 等 值 线 ,两 者 
之 间 的 斜率 与 作物 的 水 分 指数 呈 负 相关 , 据 此 提出 
了 温度 植被 干旱 指数 (TVDI) 的 概念 。TVDI 由 植被 
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图 2 TVDI 原理 示意 图 


Fig.2 Construction principle of TVDI'™) 


此 数 与 地 表 温 度 计 算得 到 ,公式 如 下 : 
T= 
TVDI = Tor (1) 
式 中 :7 为 相同 植被 指数 值 对 应 的 最 小 地 表 温 度 ， 
即 NDVI ~ T, 特征 空间 的 湿 边 ;7 为 相同 植被 指数 
对 应 的 最 大 地 表 温 度 , 即 NDVI ~ T, 特征 空间 的 干 
HT, 为 地 表 任 意 像 元 的 地 表 温 度 。 
对 最 大 最 小 地 表 温 度 散 点 进行 线性 拟 合 ,得 到 
TVDI 特征 空间 的 干 湿 边 方程 : 
Taa =a, +b, x NDVI 
Tain = @, +b, x NDVI 


式 中 ;al bi a,b, 分 别 是 干 湿 边线 性 拟 合 方程 的 系 
数 。 对 于 TVDI 值 反映 干旱 程度 的 划分 ,一 般 认 为 
值 越 趋向 于 1 ,土壤 干旱 程度 越 严重 ;TVDI 值 越 趋 
向 于 0, 土 壤 含 水 量 越 高 。 
1.3.2 ”地表 温度 的 高 程 修正 原理 ”人 研究 区 南北 高 、 
中 间 低 的 地 势 特点 ,使 得 地 形成 为 影响 研究 区 植被 、 
气候 土壤 等 地 域 分 异 的 重要 因素 。 一 方面 高 程 落 
差 较 大 的 山脉 河谷 等 地 域 格局 通过 影响 降水 .风力 
等 气候 要 素 分 布 ,对 土壤 的 含水 量 产生 影响 。 另 一 
方面 高 程 越 高 ,大 气 得 到 的 地 面 辐射 越 少 ,从 而 气温 
越 低 。 大 气 与 下 热 面 之 间 以 满 流 形式 进行 热 交 换 ， 
地 气温 差 发 生 改 变 , 通 过 热 交 换 能 量 发 生变 化 ,使 得 
地 表 温度 在 一 定 程度 上 受到 气温 的 影响 。 而 
地 表 温 度 是 TVDI 模型 的 构建 因子 ,研究 区 从 北部 
到 南部 地 形 由 高 原 过 渡 到 平原 ,又 过 渡 到 山地 , 较 大 
的 高 程 落差 会 对 地 表 温 度 产 生 一 定 的 影响 ,从 而 影 
响 TVDI 反 演 土壤 水 分 的 精度 。 

进一步 研究 TVDI 值 与 高 程 值 的 关系 ,通过 传 
统 方法 计算 陕西 省 年 平均 TVDI 值 ,并 利用 DEM 数 
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图 3 研究 区 不 同 高 程 等 级 TVDI 分 布 


Fig.3 Distribution of TVDI at different elevations in the study area 


据 提 取 每 个 高 程 值 对 应 的 TVDI 值 ,获取 研究 区 不 
同 高 程 值 对 应 的 TVDI 散 点 图 (图 3)。 从 图 3 可 以 
看 出 ,TVDI 与 高 程 呈现 明显 的 负 相 关 关系 ,相关 系 
数 为 -0.851, 说 明 研 究 区 的 TVDI 值 在 一 定 程度 上 
受到 了 高 程 的 影响 。 因 此 ,有 必要 利用 高 程 来 修正 
地 表 温 度 ,以 提高 利用 TVDI 监测 农业 干旱 的 精度 。 

本 文通 过 DEM 数据 对 地 表 温 度 进行 修正 ,所 用 
公式 如 下 : 


t, =t, +Hxm (3) 
式 中 :i 为 经 过 DEM 修正 后 的 地 表 温 度 ;i, 为 修正 
前 的 温度 ;有 厂 是 高 程 值 ;m 是 地 表 温 度 受 气温 影响 随 
高 程 增 加 而 降低 的 程度 。 
1.3.3 ”植被 指数 的 地 形 修正 原理 ”人 研究 发 现 , 同 类 
植物 处 于 阴 面 比 阳 面 同时 期 的 光谱 反射 率 低 , 尤 其 
在 近 红外 波段 差异 更 大 , 阴 面 植被 光谱 曲线 的 红 边 
斜率 随 之 降低 ,两 者 的 NDVI 值 出 现 较 大 偏差 "” 。 
根据 不 同 坡 向 的 植被 反射 率 出 现 差异 的 特征 ,相关 
学 者 认为 可 以 在 不 改变 多 光谱 曲线 特征 的 条 件 下 ， 
在 NDVI 的 计算 公式 中 采取 降低 分 母 数值 增加 分 子 
数值 来 提高 NDVI 的 值 ,同时 使 得 阳 面 比 阴 面 的 
NDVI 数值 提升 的 幅度 要 小 ,最 终 使 得 处 于 地 形 起 
伏 较 大 \ 坡 向 不 同 , 但 植被 类 型 相同 的 NDVI 值 较为 
接近 , 进而 达到 消除 或 者 抑制 地 形 影 响 的 目 
的 ”3 。 在 此 基础 上 , 吴 志 杰 等 "提出 使 用 归 一 
化 差 值 山地 植被 指数 (NDMVI) ,用 来 修正 复杂 地 形 
对 NDVI 值 的 影响 ,公式 如 下 : 
(Pair —Prea) + (Rmin — NER min ) 
(Pair Prea) — CRinin + NIR min ) 
IP pri DUETIN BE AY RE 5 prea NAL SEK E HS 
BRS FE ; R nin N ETG U Pt S OR MEL; NIR nin AU 
红外 波段 反射 率 最 小 值 。 
研究 区 高 程 落差 较 大 ,地 形 的 起 伏 变化 导致 遥 


NDMVI = 


(4) 


于 时 


感 成 像 上 出 现 不 同 的 阴影 区 ,虽然 NDVI 采用 比值 
法 可 以 抑制 部 分 由 太阳 高 度 角 不 同 造成 的 “ 同 物 异 
谱 " 现 象 , 但 地 形 的 起 伏 变 化 对 NDVI 产生 的 影响 仍 
不 能 忽视 。NDMVI 指数 是 在 近 红 外 波段 与 红 光 波 
段 分 别 减 去 一 定 的 数值 ,不 同 地 物 两 波段 间 光 谱 特 
征 的 相对 关系 并 没有 改变 ,因此 ,该 指数 可 以 使 得 阳 
面 和 阴 面 的 植被 指数 值 更 接近 ,消除 或 抑制 地 形 因 
素 造 成 的 植被 信息 偏差 。 


结果 与 分 析 


2.1 地 表 温 度 的 修正 

在 对 地 表 温 度 进行 修正 的 过 程 中 ,最 重要 的 是 
高 程 对 地 表 温 度 的 影响 程度 值 m 的 确定 。 一 般 认 
为 ,高 程 每 升 高 100 m ,气温 降低 0.6 CAA. 虽然 
随 着 气温 的 变化 ,地 表 温 度 也 会 受到 一 定 程度 的 影 
响 ,但 对 于 地 表 温 度 值 随 高 程 变化 的 研究 还 不 多 。 
相关 学 者 提出 ,相同 高 程 下 , 浅 层 地 表 温 度 值 与 气温 
值 差别 不 大 ,变化 趋势 也 相同 ,可 用 气温 随 高 程 的 变 
化 作为 参考 。 为 了 进一步 研究 高 程 对 地 表 温 度 
产生 的 影响 ,应 选取 研究 区 气温 日 较 差 相对 较 大 ,高 
程 对 气温 值 有 明显 影响 的 时 段 进 行 研究 。 研 究 区 初 
夏 气 温 日 变化 明显 ,高程 对 气温 产生 的 影响 不 容 忽 
视 , 因 此 ,选取 研究 区 6 H EARME t, 与 DEM 
数据 ,分析 地 表 温 度 随 高 程 增 加 而 降低 的 程度 mm。 
研究 区 高 程 为 171 ~3 621 m, 由 于 3 000 m 以 上 高 
程 区 域 主要 分 布 在 太白 山 等 高 山地 区 ,其 植被 类 型 
为 高 山 灌木 ,此 处 旱情 一 般 不 做 研究 ,因此 ,对 高 程 
大 于 3 000 m 的 情况 不 予 考虑 。 对 171 ~3 000 m 
的 高 程 范围 ,以 100 m 为 增 量 计算 各 个 高 程 区 间 地 
表 温 度 的 平均 值 i, 与 高 程 的 平均 值 H, ,按照 有 升 
序 排列 ,计算 相 邻 区 间 平 均 地 表 温 度 和 高 程 的 差 值 ， 
并 定义 公式 : 


2 


tay tm 
“Hp -Hn 
式 中 :mm 为 地 表 温 度 随 高 程 增 加 而 降低 的 程度 ;ti 、 
tz Hn、His 分 别 为 相 邻 区 间 地 表 温 度 的 平均 值 以 及 
高 程 的 平均 值 。 根 据 DEM 数据 提取 各 高 程 点 对 应 
的 地 表 温 度 , 并 计算 各 高 程 区 间 不 同 的 m 值 ,如 表 1 
所 示 。 
由 表 1 可 以 看 出 , m 值 波动 范围 较 小 ,大 致 分 
布 在 -0.012 71 ~0 之 间 。 根 据 获得 的 28 个 m 值 


m 


(5) 


表 1 各 高 程 区 间 的 影响 程度 值 


Tab.1 Influence values at different elevations 


高 程 区间 /m ”影响 程度 值 (m) |) 高 程 区间 /m ”影响 程度 值 (m) 
200 ~ 300 —0. 000 836 057|| 1 600 ~1 700 -0.014 751 601 
300 ~ 400 —0.001 858 027]! 1 700 ~1 800 -0.000 447 937 
400 ~ 500 -—0.003 2545 52]! 1 800 ~1 900 -0.005 699 998 
500 ~ 600 —0.001 198 907}! 1 900 ~2 000 -0.006 330 563 
600 ~ 700 —0. 007 730 485|| 2 000 ~2 100 -0.005 360 715 
700 ~ 800 —0.007 346 475| 2 100 ~2 200 -0.004 997 115 
800 ~ 900 0.002 962 541|| 2 200 ~2 300 -0.006 476 033 

900 ~1 000 -0.007 320 064]| 2 300 ~2 400 -0.004 197 357 

1 000 ~1 100 -0.010 778 05 || 2 400 ~2 500 -0.004 4604 04 

1100 ~1 200 -0.011 290 732|| 2 500 ~2 600 -0.004 246 25 

1 200 ~1 300 -0.012 257 724|| 2 600 ~2 700 -0.003 897 37 

1 300 ~1 400 -0.006 154 835|| 2 700 ~2 800 -0.009 9739 66 

1 400 ~1 500 -0.012 229 793 2 800 ~2 900 -0.002 819 414 

1 500 ~1 600 -0.012 709 305|| 2 900 ~3 000 -0.002 372 969 


分 别 对 6 月 上 旬 的 地 表 温 度 进行 修正 ,利用 修正 后 
的 地 表 温度 数据 结合 NDVI 值 反 演 TVDI, 并 与 分 布 
于 人 研究 区 28 个 农业 气象 站 点 测 得 的 相应 时 间 的 20 
cm 土壤 相对 湿度 数据 做 相关 性 分 析 。 如 反 演 得 到 
TVDI 值 与 实测 土壤 相对 湿度 数据 有 最 好 的 相关 性 ， 
则 反 演 得 到 的 TVDI 值 能 最 好 的 反映 土壤 湿度 状 
况 ,说 明 利用 该 m 值 修正 后 的 地 表 温度 能 够 较 准 确 
地 反映 地 表 实际 温度 , 则 取 为 最 佳 m 值 。 如 图 4 所 
示 , 由 于 气温 与 高 程 呈 现 明显 的 负 相 关 关 系 ,为 了 能 
与 气温 随 高 程 的 变化 相对 应 ,并 使 得 计算 表达 方便 ， 
Ki m RIRN -1/100 为 单位 的 形式 。 将 利用 不 同 m 
值 修正 后 的 地 表 温度 有 反 演 得 到 的 TVDI 值 与 土壤 湿 
度数 据 进行 相关 性 分 析 ,其 分 布 如 图 4。 

根据 影响 程度 m 与 相关 系数 r 的 关系 图 可 知 ， 


1.00 p 
0.90 + 
te 0.80 | 
Ex 0.70 + 
BRO 60 | 
HR 
ip tK 0.50 p 
IE 0.40 + 
Si 0.30 + 
fa BA 0.20 | 
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Fig.4 Relationship between the influence level of elevation 


图 4 


on surface temperature and the correlation coefficient 
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随 着 m 值 的 逐渐 增 大 , 反 演 得 到 的 TVDI 值 与 实测 
土壤 相对 湿度 的 相关 系数 呈现 先 增加 后 减少 的 趋 
势 , 晶 最 适合 的 m 值 落 在 -0. 45/100 ~ -0. 75/100 
之 间 。 进 一 步 分 析 最 佳 m 值 ,发 现在 -0. 45/100 ~ 
-0. 75/100 区 间 相 关系 数 r 呈现 明显 的 先 增加 后 降 
低 的 趋势 ,有 目 在 m= -0.006 15 时 相关 系数 达到 最 
大 ,此 时 反 演 得 到 的 TVDI 与 实测 土壤 湿度 数据 有 
最 好 的 相关 性 , 且 这 个 值 与 气温 随 着 高 程 增 加 而 降 
低 的 程度 相近 ,说 明 取 影响 程度 m 为 -0.006 15 是 
可 信 的 ,由 此 得 到 高 程 对 地 表 温 度 的 影响 程度 最 佳 
(E m= -0.006 15。 根 据 求 得 的 最 佳 m 值 ,对 研究 
区 2010 年 4 月 7 日 至 9 月 30 日 地 表 温 度 分 别 进行 
修正 ,以 减弱 高 程 对 地 表 温 度 的 影响 。 
2.2 植被 指数 的 修正 

为 了 减弱 地 形 起 伏 对 获取 植被 信息 的 影响 ,本 
研究 采用 NDMVI 提取 植被 信息 。 根 据 植 被 指数 计 
算 公 式 分 别 求 得 陕西 省 2010 年 4 月 7 日 至 9 月 30 
H 16 d 合成 的 NDVI 与 NDMVI 影像 ,为 了 更 加 明显 
地 对 NDVI 与 NDMVI 影像 进行 对 比 评价 ,在 保证 数 
据 源 相同 的 基础 上 ,选择 太阳 光照 强度 适中 的 4 月 
23 HÉJ NDVI 与 NDMVI 影像 进行 分 析 ( 图 5)。 从 
图 5 可 以 看 出 ,NDVI 与 NDMVI 影像 植被 指数 值 总 
体 分 布 相似 ,主要 区 别 在 于 由 地 形 起 伏 变化 形成 的 
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阴影 区 与 非 阴影 区 。 从 图 Sa 可 以 得 出 :影像 表现 为 
明显 的 地 形 纹理 特征 , 且 在 地 形 起 伏 变 化 较 大 的 阴 
影 区 域 影像 NDVI 值 较 低 。 而 图 Sb 可 以 看 出 ,在 
NDVI 影 像 上 植被 指数 值 较 低 的 区 域 NDMVI 值 有 
所 增加 , 且 整 体 影像 地 形 纹理 特征 不 明显 ,在 地 形 起 
伏 变 化 较 大 的 区 域 植被 指数 值 有 所 提高 。 可 见 ， 
NDMVI 使 得 由 地 形 起 伏 变化 造成 的 阴影 区 与 非 阴 
影 区 的 植被 信息 较为 接近 ,可 以 帮助 消除 或 者 抑制 
由 地 形 差异 而 导致 植被 信息 失真 的 现象 。 
2.3 TVDI 指数 计算 结果 

本 文 利用 经 过 高 程 和 地 形 修正 前 后 的 地 表 温 度 
和 植被 指数 分 别 构建 4 种 TVDI 指数 :QD 利用 原始 
地 表 温 度 与 NDVI 构建 的 TVDI, ;@@ 利用 经 过 DEM 
修正 后 的 地 表 温 度 与 NDVI WER TVD aO 利用 
原始 地 表 温 度 与 NDMVI 构建 的 TVDI,,; 由 利用 
DEM 修正 后 的 地 表 温 度 与 NDMVI 构建 的 TVDL,。 

将 植被 指数 值 进 行 统 计 分 析 , 以 0.01 间隔 为 步 
长 ,提取 相同 植被 指数 值 对 应 的 最 大 最 小 地 表 温 度 。 
根据 最 大 、 最 小 地 表 温 度 的 散 点 图 进行 线性 拟 合 , 即 
得 到 干 湿 边 方程 。 结 果 见 表 2 ~5。 

从 TVDL 与 TVD1,, 的 干 湿 边 拟 合 方程 结果 来 
看 ,在 12 期 影像 中 有 10 期 TVDI,, 的 干 边 方程 RR 较 
TVDI, 有 所 提高 ,其 均值 由 0.786 提 高 至 0.817 ; 湿 


(b) NDMVI 


5 NDVI/NDMVI 指数 影像 
Fig.5 Images of NDVI(a) and NDMVI(b) 


边 方程 RR 有 6 期 影像 较 TVDL, 的 湿 边 方程 R 得 到 
了 提高 ,其 均值 由 0. 219 提高 至 0.305。 由 此 得 知 ， 
TVDI,, 的 干 湿 边 散 点 分 布 较为 聚集 ,而 TVD 的 干 
湿 边 散 点 分 布 较为 离散 ,说 明 TVDI, 的 干 湿 边 方程 
拟 合 效果 较 好 ,经 过 高 程 修正 后 的 地 表 温 度 更 能 准 
确 反 映 实 际 地 表 温 度 的 分 布 状况 。 从 TVDI, 与 


2.4 4 种 TVDI 指数 的 比较 与 验证 

2.4.1 TVDI 与 土壤 湿度 的 关系 ”为 了 对 比分 析 4 
种 TVDI 监测 农业 旱情 的 效果 ,笔者 计算 了 陕西 省 
2010 年 4 一 9 月 4 种 TVDI 值 ,同时 提取 各 土壤 湿度 
站 点 对 应 的 4 种 遥感 指数 值 ,并 结合 相应 时 间 的 实 
测 土壤 相对 湿度 数据 进行 相关 性 分 析 ( 图 6)。 结 


TVDI, 的 干 湿 边 拟 合 方程 结果 来 看 ,TVDI, 的 拟 合 
效果 整体 较 好 ,其 中 干 边 方程 RR 均 较 TVDI, 有 所 提 
高 ,其 均值 提高 至 0. 891; 湿 边 方程 R 有 8 期 影像 
都 有 所 提高 ,其 均值 由 0.219 提高 至 0.331, 说 明 
NDMVI 指数 能 更 加 准确 地 获取 地 表 植 被 信息 。 同 
时 ,对 比分 析 TVDL 与 TVDI ,的 干 湿 边 拟 合 方程 可 
知 ,TVDIu 的 干 湿 边 方程 拟 合 结果 较 好 ,其 中 干 边 方 
FE R 均值 由 0.786 提高 至 0.915 , 湿 边 方程 尼 均值 
由 0.219 提高 至 0.339 ,说 明 TVDI ,指数 能 更 加 准 
确 的 反映 地 表 真 实 信息 。 


显示 ,4 种 温度 植被 干旱 指数 值 均 与 实测 土壤 相对 
湿度 数据 呈现 负 相关 关系 。 

由 图 6 分 析 可 知 ,TVDI,, 值 与 土壤 相对 湿度 的 
相关 性 高 于 TVDL, 值 ,其 相关 系数 > 提高 了 0. 204, 
说 明 经 过 高 程 修正 后 的 TVDIL, 能 较 好 地 反映 真实 
土壤 湿度 的 状况 ,可 提高 干旱 监测 的 精度 。 通 过 对 
比 TVDI, 和 TVDI, 值 与 土壤 湿度 的 关系 可 知 ,两 者 
均 与 土壤 湿度 呈现 负 相 关 关 系 ,但 TVDL, 值 与 土壤 
相对 湿度 相关 性 更 好 ,相关 系数 为 -0.460, 较 原始 
TVDI, 值 与 土壤 相对 湿度 的 相关 系数 提高 了 0. 118。 


表 2 TVDL, 指数 的 干 湿 边 方程 
Tab.2 The dry-wet edge fitting equations for TDVI, 


日 期 /月 -日 干 边 拟 合 方程 湿 边 拟 合 方程 R? 
04 -07 y = -12.005x + 27.633 0. 826 9 y = 7.1061x - 1.4405 0. 158 7 
04 -23 y = -32.92x + 37.615 0.920 1 y = 19. 546x - 5.174 0.735 8 
05 -09 y = —35. 464x + 42.825 0. 888 4 y = 6.7253x — 1.986 1 0. 304 6 
05 -25 y = -38.757x + 43.93 0. 856 3 y = 2.331 7x + 0.0196 0.021 7 
06 - 10 y = -23.067x + 45.117 0. 806 2 y = -13.111x + 10.148 0.295 6 
06 -26 y = —25.086x + 49.516 0.757 3 y = —8.481 8x + 16.522 0.180 9 
07 -12 y = -20.612x + 44. 853 0.749 9 y = 0.651 7x + 0.908 0.008 1 
07 -28 y = -31.384x + 46.476 0. 708 7 y = 5.203 9x - 1.4696 0. 234 0 
08 -13 y = —18.778x + 36. 628 0. 734 8 y = -16.025x + 8.9533 0.410 1 
08—29 y = —24.896x + 36.462 0. 825 0 y = -8.458 lx + 6.3947 0. 076 9 
09 - 14 y = —30.486x + 38.177 0. 602 0 y = 8.896 8x - 2.6735 0. 206 6 
09 -30 y = -32.281x + 36.289 0. 752 6 y = -0.707 7x + 1.0187 0. 003 2 

#3 TVDL,, 指 数 的 干 湿 边 方程 
Tab.3 The dry-wet edge fitting equations for TVDI,, 

日 期 /月 -日 干 边 拟 合 方程 湿 边 拟 合 方程 R? 
04 -07 y = -14.505x + 29.074 0. 883 1 y = 11. 054x - 3.7167 0. 246 2 
04 -23 y = -21.318x + 36.063 0. 962 0 y = -2.696 9x + 12.645 0. 209 2 
05 -09 y = —24.254x + 41.877 0. 848 6 y = -14.479x + 11.476 0.315 7 
05 -25 y = -21.942x + 40.792 0. 900 6 y = -4.550 5x + 4.3192 0.091 9 
06 -10 y = —28.617x + 48.904 0. 927 6 y = -5.129 4x + 4.6845 0.175 4 
06 -26 y = —25.929x + 49.396 0.765 2 y = -7.051 5x + 16.163 0.1100 
07 - 12 y = —20. 834x + 44.696 0. 769 8 y = -0.021 7x + 1.3343 0. 436 1 
07 -28 y = —16.787x + 42. 858 0.901 7 y = 5.026 9x — 1.627 0. 188 2 
08 - 13 y = —14.294% + 36.651 0. 707 7 y = —18.64x% + 16.609 0. 664 8 
08 -29 y = -14.017x + 34.374 0. 684 4 y = -18.417x + 16.284 0.559 5 
09 - 14 y = -15.781x + 34.565 0. 693 1 y = -1.308 5x + 1.8781 0.039 1 
09 -30 y = -17.764x + 33.202 0.755 9 y = -13.026x + 10.221 0. 620 9 
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表 4 TVD, 指数 的 干 湿 边 方程 
Tab.4 The dry-wet edge fitting equationsfor TVDI,, 


日 期 /月 -日 干 边 拟 合 方程 R? 湿 边 拟 合 方程 R? 
04 -07 y = -14. 684% + 33.224 0. 871 8 y = 14.416x + 0.623 0. 600 5 
04 -23 y = -28.428x + 38.873 0. 963 8 y = 24.302x + 3.4809 0. 824 1 
05 -09 y = -34.967x + 46.518 0. 896 6 y = 0.725 1x + 12.235 0. 003 6 
05 -25 y = -37.92x + 48.683 0.969 5 y = 7.888 4x + 8.298 3 0.291 0 
06 - 10 y = -22.325x + 50.049 0. 865 3 y = -5.890 9x + 17.29 0. 170 6 
06 -26 y = -24.009x + 53.395 0. 827 0 y = -6.194 7x + 24.589 0.185 5 
07 -12 y = -25. 184x + 53. 867 0. 874 8 y = 2.333 1x + 12.48 0. 088 0 
07 -28 y = -38.618x + 55.400 0.8713 y = 10. 871x + 8.7507 0.574 8 
08 -13 y = -26.605x + 44. 126 0.911 0 y = -7.536 4x + 18.859 0.344 9 
08 -29 y = -30. 14x + 43.815 0.963 4 y = 3.495 8x4 11.501 0.027 3 
09 -14 y = -23.996x + 38. 862 0.793 8 y = 15.382x + 4.3319 0. 454 9 
09 -30 y = -25.423x + 36.434 0. 887 5 y = 14.409x + 2. 188 5 0. 416 8 

#5 TVDI,, 指 数 的 干 湿 边 方程 
Tab.5 The dry-wet edge fitting equations for TVDI,,, 

日 期 /月 -H 干 边 拟 合 方程 R? 湿 边 拟 合 方程 R? 
04 -07 y = -17.363x + 34.764 0. 930 6 y = 17.978x - 1.4341 0. 652 7 
04 -23 y = -19.868x + 38.725 0. 985 8 y = 13.241x + 4.2838 0.710 8 
05 -09 y = -22.921x + 45.001 0. 945 3 y = -5.988 3x + 17.126 0.321 2 
05 -25 y = -26.009x + 48.446 0. 981 0 y = -1.050 6x + 13.126 0.075 0 
06 -10 y = -26.896x + 53.173 0. 950 0 y = -2.686x + 15.101 0. 263 2 
06 -26 y = -25.633x + 53.623 0. 856 2 y = -3.158 9x + 23.082 0.345 0 
07 -12 y = -26.353x + 53.760 0.923 6 y= 1.960 6x +12.974 0. 064 8 
07 -28 y = -23.614x + 51.714 0.967 1 y = 7.246 8x + 9.858 2 0.398 2 
08 -13 y = -17.427x + 42.705 0. 894 6 y = -9.064x + 22.687 0. 568 1 
08 -29 y = -18.075« + 41.971 0. 884 2 y = -10. 105x + 20. 805 0. 336 0 
09 -14 y = -15.131x + 37.437 0.770 6 y = 1.125 6x + 10.449 0. 024 0 
09 -30 y = -15.541x + 34. 993 0. 887 1 y = -8.926 5x + 13.508 0. 314 3 

120 p (a) NDVI-T, r= -0.342 120 p (b) NDVI-T (DEM) r=- 0.546 
& 100 上 os P<001 Pees (|e p P<0.01 
其 oe & a Peet ee ‘en -3 
m 80 f R ge 80 + * 
a 60 上 E 60 + 
I 40 上 Æ 40 + 
H 20| E 20 上 
0 i 1 4 0 | j 
0.20 0.40 0.60 0.80 1.00 0.20 0.40 0.60 0.80 1.00 
TVDI, TVDI,, 
, 120; ()NDMVI-T, r= -0.460 120 | (d) NDMVI-T,(DEM) r=-0.711 
a 100 上 x= 100 + 
mj 80 上 由 80 + 
me 60 | = 60 | 
40 上 = 40} 
H 20} H 20 | 
0 J 0 f 
0.20 0.40 0.60 0.80 1.00 0.20 0.40 0.60 0.80 1.00 
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图 6 4 种 TVDI 与 土壤 相对 湿度 的 关系 
Fig.6 Relationship between the four kinds of VDI and the soil relative humidity 
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图 7 降水 量 与 TVDI 在 站 点 尺度 随时 间 变 化 关系 


Fig.7 Relationship between precipitation and TVDI on station scale 
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郑 发 美 等 :基于 地 形 因子 的 TVDI 修正 一 一 以 陕西 省 为 例 


从 修正 前 后 的 植被 指数 值 与 土壤 相对 湿度 的 相关 性 


可 以 减弱 因 高程 引 起 的 相同 植被 覆盖 条 件 下 ,高程 


来 看 ,NDMVI 可 以 更 加 准确 地 反映 地 表 植 被 信息 ， 
对 比 原始 TVDI, 和 修正 后 的 TVDL,, 值 与 土壤 相对 
湿度 的 相关 性 可 知 ,TVDL , 值 与 土壤 相对 湿度 的 相 
关 性 更 好 ,相关 系数 从 - 0. 342 提高 至 -0.711。 修 
正 后 获取 的 TVDL. 值 随 实测 土壤 相对 湿度 的 增加 ， 
呈现 明显 降低 的 趋势 ,说 明 TVDI, 能 更 好 地 反映 土 
壤 湿 度 状况 ,能 更 加 准确 地 进行 干旱 监测 。 

2.4.2 TVDI 与 降水 的 关系 ， 在 气象 站 点 尺度 上 ， 
探讨 4 种 TVDI 值 与 降水 的 关系 。 按 照 空间 分 布 均 
匀 性 .站 点 所 在 像 元 成 像 清晰 性 等 原则 ,选择 榆林 、 
横山 延安 等 12 个 气象 站 点 数据 进行 分 析 。 为 研究 
4 种 TVDI 值 与 降水 的 同步 性 ,计算 各 个 站 点 对 应 各 
期 遥感 影像 时 间 内 的 月 降水 均值 ,同时 提取 各 期 影 
像 中 每 个 气象 站 点 处 4 种 月 尺度 的 TVDI 均值 。 分 
析 各 TVDI 值 与 降水 的 变化 关系 (图 7) 。 

从 图 7 可 以 看 出 ,在 4 一 9 月 时 间 段 内 ,4 种 TV- 
DI 值 的 变化 趋势 总 体 上 一 致 , 且 波 动 均 与 降水 波动 
特征 相反 。 在 波动 幅度 上 ,原始 TVDI 较 修 正 后 的 3 
种 TVDI 数值 波动 幅度 较 弱 ,说 明 经 过 修正 的 TVDI 
较 原始 TVDI 对 降水 的 响应 更 强 , 即 其 对 旱情 的 变 
化 更 加 敏感 。 而 在 降水 量变 化 较 大 的 8 月 ,TVDI 值 
在 陕 北 SEP DERG 3 个 地 区 呈现 出 不 同 的 波动 特 
征 。 其 中 , 陕 北 地 区 的 榆林 、 定 边 .横山 ` 绥 德 及 延安 
5 个 降水 记录 站 点 均 体现 出 VDI, > TVDI,, > 
TVDI, > TVDL 的 特征 ,说 明 陕 北 地 区 地 形 起 伏 较 
大 ,经 过 高 程 修 正 后 的 TVDI,, 与 TVDI,, 值 更 能 准确 
反映 土壤 水 分 售 量 ;而 在 关中 地 区 的 长 武 . 洛 川 , 滩 
县 及 武功 4 个 降水 记录 站 点 均 呈 现 出 TVD aa > 
TVDL,, > TVDIL > TVDL, 的 特征 ,说明 关 中 地 区 地 势 
相对 平坦 ,但 由 修正 后 的 植被 指数 计算 得 到 的 
TVDL,, 与 TVDI,, 值 与 降水 的 波动 更 为 贴切 ; 陕 南 地 
DOA FE I AR 3 个 降水 记录 站 点 体现 出 
TVDI,,, > TVDI,, > TVDI, > TVDI 的 特征 , 即 修正 后 
的 TVDI 对 月 尺度 降水 的 响应 要 好 于 未 做 修正 的 
TVDI, H. TVDIL. 最 优 。 绽 上 ,经 过 高 程 和 地 形 起 伏 
修正 后 获取 的 TVDI,, 能 较 好 的 反映 出 降水 的 变化 
特征 ,对 旱情 有 较 好 的 监测 效果 。 


3 结论 


(1) 高 程 通过 对 地 表 温 度 的 影响 ,进一步 影响 
TVDI, 利 用 DEM 数据 对 地 表 温 度 进行 了 地 形 校正 ， 


高 的 像 元 对 应 的 地 表 温 度 值 降低 、TVDI 减 小 的 现 
Zo Ait DEM 数据 修正 后 的 TVDI,, 值 与 土壤 湿度 
的 相关 系数 较 原 始 TVDI, 提高 了 0. 204, 说 明 修 正 
后 的 地 表 温 度 更 能 准确 反映 实际 情况 。 

(2) 与 NDVI 相 比 ,NDMVI 能 够 有 效 降低 因 地 
形 起 伏 对 植被 信息 提取 产生 的 影响 。 经 过 地 形 修正 
后 的 TVDI, 值 与 土壤 相对 湿度 的 相关 系数 较 原 始 
TVDI, 值 提 高 了 0. 118。 表 明 NDMVI 可 以 更 加 准确 
地 反映 植被 信息 ,在 一 定 程度 上 可 以 消除 地 形 影 响 
和 削弱 地 形 阴影 的 干扰 。 

(3) 经 过 高 程 和 植被 指数 修正 后 的 TVDI,。 与 
土壤 相对 湿度 数据 的 相关 系数 为 -0.711, 较 原始 
TVDIn 提高 了 0. 369 ,说 明 修 正 后 的 TVDI,, 能 更 好 
地 反映 土壤 湿度 状况 。 同 时 ,TVDT, 与 降水 波动 特 
征 具 有 较 好 的 对 应 关系 ,说 明 经 过 地 形 因 子 修正 后 
的 TVDI 可 对 农业 旱情 有 较 好 的 监测 效果 。TVDL,。 
可 有 效 减弱 高 程 与 地 形 起 伏 对 TVDI 指数 的 影响 ， 
且 原 理 简 单 ,计算 方便 ,效果 较 TVDI, 显著 提高 ,不 
受 其 他 条 件 限 制 ,具有 较 高 的 实用 意义 。 
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Abstract: In recent years , Temperature Vegetation Drought Index (TVDI) has been widely applied in monitoring 
drought with remote sensing. However, the influence of topographical factors on the index accuracy is not considered 
in the traditional TVDI. In this paper,a case study on the TVDI modification in Shaanxi Province was carried out so 
as to reduce the influence of topographical factors on the monitoring accuracy of TVDI,the DEM data were used to 
make an elevation correction of surface temperature ,the NDVI was replaced by the Normalized Difference Mountain 
Vegetation Index (NDMVI) so as to reduce the influence of topographical fluctuation on vegetation index, and the 
new TVDI indexes were constructed. Four kinds of TVDI indices were constructed by using the surface temperature 
and NDMVI and NDVI before and after correction, and the relationships among the 4 indexes and soil moisture and 
precipitation were analyzed so as to compare the accuracy of monitoring agriculture drought. The results showed that 
the elevation and terrain fluctuation affected significantly the TVDI precision in monitoring drought, and the TVDI 
after the correction of surface temperature and vegetation was better for monitoring drought compared with the tradi- 
tional TVDI, its correlation coefficient with relative soil humidity was increased from 0.342 to 0.711, and it could 
be used to accurately reflect the change of precipitation. The results revealed that the modified TVDI could be used 
to effectively reflect the soil water conditions ,and could also be used to accurately monitor agricultural drought. 
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